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Abstract 

We present the latest Next-to-Next-Leading order strong coupling constant, Og, ex- 
tracted from Gross-Llewellyn-Smith sum rule. The value of Og from this measurement, at 
the mass of Z boson, is a^^'^'^(M|) = 0.114to;o?2- We discuss the previous strange sea 
quark measurement from the CCFR experiment and the prospects for improvements of 
the measurement in the NuTeV experiment. 



Introduction 



Neutrino-nucleon (u-N) deep inelastic scattering (DIS) experiments provide a good test- 
ing field for Quantum Chromo-Dynamics (QCD), the theory of strong interactions. The u-N 
DIS experiments probe the structure of nucleons and provide an opportunity to test QCD 
evolutions and to extract the strong coupling constant, a.^- They are complementary measure- 
ments to charged lepton DIS experiments of nucleon structure functions. The advantage of u-N 
DIS measurements over charged lepton experiments is that z/-N experiments can measure both 
nucleon structure functions, F2{x,Q'^) and xF^lxjQ^), due to pure V-A nature. The u — N 
differential cross sections are written, in terms of the structure functions: 

\ ^ 2 l + R{x,Q^) I ^ V 2 / J 

where R{x, Q"^) = (Jl/(^t-, the ratio of longitudinal to transverse absorption cross sections. The 
structure functions F2{x, Q^) and xF^{x, Q"^) are extracted from fitting Eq. [l| to the measured 
differential cross sections. 

In this paper, we present the next-to-next leading order (NNLO) determination of as 
from the Gross-Llewellyn-Smith (GLS) sum rule which states that the total number of valence 
quarks are given by the integration of the non-singlet structure function xF^, over entire regions 
of x, the momentum fraction carried by the struck quark. We also discuss the measurements 
of strange sea quark distributions from CCFR experiment and the expected improvements of 
this measurement in NuTeV. 

The Experiment 

CCFR (Columbia-Chicago- Fermilab-Rochester)/NuTeV experiment is a. v — N DIS ex- 
periment at the Tevatron in Fermilab. The CCFR experiment used broad momentum beam 
of mixed neutrinos (z/^) and antineutrinos iv^) from decays of the secondary pions and kaons, 
resulting from interactions of 800 GeV primary protons on a Beryllium- Oxide (BeO) target. 
The NuTeV experiment, successor of CCFR, used separate beams of z/^ or during a given 
running period, using a Sign-Selected-Quadrupole- Train (SSQT) [0. 

The CCFR/NuTeV detector § consists of two major components : target calorimeter 
and muon spectrometer. The target calorimeter is an iron-liquid-scintillator sampling calor- 
imeter, interspersed with drift chambers to provide track information of the muons result- 
ing from charged- current (CC) interactions where a charged weak boson (W"*" or W~) is ex- 
changed between the i/^ (l7^) and the parton. The calorimeter provides dense material in 
the path of v iv) to increase the rate of neutrino interactions. The hadron energy reso- 
lution of the calorimeter is measured from the test beam and is found to be : a / Euad = 

(0.847 ± Qmb)/^EHad{GeV) + (0.30 ± ().12) /Euad- 

The muon spectrometer is located immediately downstream of the target calorimeter 
and consists of three toroidal magnets and five drift chamber stations to provide accurate 
measurements of muon momenta. The momentum resolution of the spectrometer {a/p^) is 
approximately 10.1% and the angular resolution is 6^ = 0.3 + 60/p^(mrad). 

as from Gross-Llewellyn- Smith Sum Rule 

Once the structure function xF^ is extracted, one can use Gross-Llewellyn-Smith (GLS) 
sum rule 0, which states that / xF^{dx/x) is the total number of valence quarks in a nucleon, 
up to QCD corrections, to extract the strong coupling constant, a^. Since in leading order 
(LO), the structure function xF^ is xq — xq, the valence quark distributions, integrating xF^ 
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Figure 1: xF^ vs x for four lowest bins. The solid circles represent the CCFR xF^ data and 
the inverse triangle represent the data from other experimental measurements. 

Since the GLS sum rule is a fundamental prediction of QCD and the integral only 
depends on valence quark distributions, as can be determined without being affected by less 
well known gluon distributions. Moreover, since there are sufficient number of measurements of 
xF^ in a wide range of Q^, one can measure as a function of in low where the values 
of as depends most on Aqcd- 

With next-to-next-to-leading order QCD corrections Q, the GLS integral takes the form 

/' xF,{x, Q^)- = 3 (l - ^ - a{n,){^r - b{nf){^f) - AHT (2) 

where the term AHT is the corrections from higher-twist effects. The higher-twist correction 
term, AHT, is predicted to be significant in some models f^, while others predict negligibly 
small corrections. We take AHT as one half the largest model prediction with the associated 
error covering the full range {AHT = (0.15 ± 0.15)/Q^). 

In order to perform the integration of xF^ in the entire ranges of x, we use data from 
other I/-N DIS experiments (WA59, WA25, SKAT, FNAL-E180, and BEBC-Gargamelle §), to 
cover large-x ranges that are not covered by CCFR due to geometric and kinematic acceptances. 
The CCFR data has a minimum x of roughly x = 0.002Q^. To extrapolate below this kinematic 
limit, we fit xF^ to a power law (Ax^) using all the data points in x < 0.1. 

Since the data points do not populate the region x > 0.5 finely, it is necessary to 
interpolate between the data points. Thus, we use the shape of the charged lepton DIS F2 in 
X > 0.5, because the shapes of xF^ and F2 should be the same in this region of x due to negligible 
contribution from sea quarks. The charged lepton F2 data are corrected for nuclear effects 
before obtaining the shapes. Systematic uncertainties due to extrapolation and interpolation 
take into account the possible variations of the models in the low-x region and the resonance 



Figure |1| shows xF^ in the four lowest bins as a function of x. The sohd circles 
represent the experimental xF^ data, The solid lines represent the power law fit used for ex- 
trapolation of X-F3 outside the kinematic limit of CCFR. The shaded area in each plot shows 
the region of x where the extrapolations are used for the integration. 

The xF^ integrals are estimated in six bins. The results in each bin are fit to 
the NNLO pQCD function and higher-twist effect in Eq. |^. This procedure yields a best fit 
of K^^^^^ = 165MeV. Evolving this result to the mass of the Z boson, Mz, in NNLO, this 
corresponds to the value of a^: 

a^NLO^Ml) = 0.114+°:°°^(.tat.);no9(^2/^i-) ± 0.005(tWy) (3) 

where systematic uncertainty includes the error in the fit and theory error represents the uncer- 
tainty due to higher twist effects and higher order QCD contributions. This result is consistent 
with the value extracted from the CCFR structure function measurement 

Evolving to the mean value of (= 3GeV^) for this analysis results in Q!^^'^*^(3GeV^) = 
0.278±0.035ib0.05l[]:|jf . If the higher-twist effect is neglected, the value of ^■^°(M|) becomes 
0.118. 

Strange Sea Measurements 

In u — N DIS, the flavor changing weak CC interactions provide a clean signature for 
scattering off an s quark (s c), resulting in a pair of oppositely signed muons in the final 
state through the reaction I'fj.iVfi) + N ^ +c(c) +X where c(c) subsequantly undergoes 

a semileptonic decay c(c) s{s) + ji^^ji^) +V^{u^). In contrast, detecting s quark with NC 
interactions in a charged lepton DIS is convoluted with s production and fragmentation of 
strange mesons, requiring good particle identifications. The nucleon structure function F2 from 
V — N DIS together with that measured from charged lepton DIS could also measure s quark 
distributions. In addition, the difference of the v — N structure function xF^ between v and V 
measures the nucleon strange quark contents as well. 

CCFR experiment performed both LO |Ty] and NLO analyses |jTl[]. Since a neutrino 



from the semileptonic decay of the charm is involved in the final state, the CCFR measured 

2 

the distributions of visible physical quantities; Xms, Ems = Enad + pi +pL Zms = p 2 
which is the fractional momentum taken by the charm meson fragmentation. The measured 
distributions were compared to Monte Carlo predictions which incorporate detector acceptance 
and resolutions. 

The Monte Carlo included the singlet {xqsi{x, /i^)), non-singlet {xq^s^x, ij?) = xqv{x, /i^)), 
and gluon parton distribution functions from CCFR structure function analysis assuming: 

xqv{x, /i^) = xuvix, fi^) + xdvix, yU^) (4) 

xdv{x, /x^) = Aci{l — x)xuvix, f/) (5) 

xu{x, jj?) = xus{x^ H^) (6) 

xd{x, fi"^) = xds{x, fi"^) (7) 

The above assumes symmetric sea quark distributions. We then parameterized the strange sea 
distribution to: 

xu{x, /i^) + xd{x, jj?) 



XS{X, jJL ) = As{l — x)° 



(8) 



where, Ag is defined in terms of level and shape parameters, n and a. The level parameter n is 
defined as : 

_ Jq[xs{x,ii^) +xs{x,ii^)] 



Table 1: Summary of the CCFR strange sea analysis. The first set of the errors in each item 
is statistical and the second set is the systematic uncertainties. 





i'ragmentation 




a 


Be 


mc(GeVj 




Collins-Spiller 


0.477 


-0.02 


0.1091 


1.70 


NLO 


e = 0.81 
±0.14 


+0.046 
-0.044 
+0.023 
-0.024 


+0.60 
-0.54 
+0.28 
-0.26 


+0.0082 
-0.0074 
+0.0063 
-0.0051 


±0.17 

+0.09 
-0.08 




Peterson 


0.468 


-0.05 


0.1047 


1.69 


NLO 


e = 0.20 


+0.061 
-0.046 


+0.46 
-0.47 


±0.0076 


±0.16 




±0.04 


+0.024 
-0.025 


+0.28 
-0.26 


+0.0065 
-0.0052 


+0.12 
-0.10 




Peterson 


0.373 


2.50 


0.1050 


1.31 


LO 


e = 0.20 


+0.048 
-0.041 


+0.60 

-0.55 


±0.007 


+0.20 
-0.22 




±0.04 


±0.018 


+0.36 
-0.25 


±0.005 


+0.12 
-0.11 



The LO prediction included a simple parton model which has tree level calculations for 
scattering of a off an s quark, resulting in a charm quark final state. On the other hand, 
the NLO prediction included the ACOT formalism [0 whose calculation is performed in a 
variable flavor MS scheme. Specifically the NLO prediction included the Born and gluon-fusion 
diagrams. In addition, since the charm quark is massive, the threshold effect was taken into 
account in the prediction via leading order slow rescaling formalism where the scaling variable 
X is replaced with ^ defined as : 

e = Ml + 7^1-^K (10) 



where ttt-c is the mass of the charm quark. 

We then fit the Monte Carlo to data distributions of x^s^ E^is, and z^s for k, a, Be (the 
charmed meson semi-leptonic branching ratio), rric, and e (fragmentation parameter), using the 



values of CKM martix elements | Vd \ and | V^s |, taken from PDG [|13|. 

Table |1| summarizes the results of this analysis, comparing various fragmentation func- 
tions at LO and NLO. The result shows that the strange sea level parameter, k, is in qualitative 
agreement between LO and NLO analyses while the shape parameter, a, is the same as other 
sea quarks in NLO but softer in LO. The charm quark mass differs between NLO and LO, 
due presumably to the fact that the parameter rric is not a pole mass but rather a theoretical 
parameter that absorbs the lack of higher corrections. 

Aside from the large statistical uncertainty, there were two major sources of systematic 
uncertainties in performing this measurement in CCFR. First and foremost one is the identifi- 
cation of the secondary muons resulting from the decay of charmed mesons. Since the neutrino 
beam in CCFR was a mixture of z/^ and V^, the experiment relied on transverse momentum of 
the muons relative to the neutrino beam direction, PJ^, to distinguish the prompt muons from 
secondary muons. The second systematic uncertainty is the LO slow rescaling formalism to 
take into account the charm threshold effect in theoretical predictions. 

The experimental systematic uncertainty due to the identification of the secondary muon 
almost completely disappears in the NuTeV experiment due to the use of SSQT whose wrong 
sign neutrino contamination at a given mode is less than 10"'^. In addition, the NuTeV ex- 
periment is working on improving low energy muon energy measurements in the calorimeter 



ning to incorporate more complete NLO calculations and will investigate the dependence of the 
predictions on factorization scheme and various parton distribution function parameterizations. 

While statistical uncertainty dominated in the previous measurement, improving the 
systematic uncertainty would reduce the total error significantly. We expect the statistical 
uncertainty in this analysis would reduce by about a factor of 2.2 combining CCFR and NuTeV 
data. The NuTeV is currently planning to present a LO analysis this summer. 

Conclusions 

CCFR measured GLS integral in six bins and extracted the value of NNLO at 
the mass of the Z boson. The measured value of is: 

«™^(M|) = (11) 

which is consistent with the value measured from the structure function analysis and with 
world average. 

NuTeV has finished data taking in September 1997, running with separate beams of 
and 17^, and a precision calibration (0.3%) is underway. Strange sea measurement at the 
NuTeV experiment is making progress, taking advantage of the separate beam and the intensive 
calibration. NuTeV expects a LO result in the strange sea measurement by the summer. 
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